Electronic structure of overstretched DNA 
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Minuscule molecular forces can transform DNA into a structure that is elongated by more than half 
its original length. We demonstrate that this pronounced conformational transition is of relevance 
to ongoing experimental and theoretical efforts to characterize the conducting properties of DNA 
wires. We present quantum mechanical calculations for acidic, dry, poly (CG) -poly (CG) DNA which 
has undergone elongation of up to 90 % relative to its natural length, along with a method for 
visualizing the effects of stretching on the electronic eigenstates. We find that overstretching leads to 
a drastic drop of the hopping matrix elements between localized occupied electronic states suggesting 
a dramatic decrease in the conductivity through holes. 

PACS numbers: 87.14.Gg, 72.80.Le, 87.15.Aa 



The possibility of efficient charge transfer along the 
stacked tt orbitals in DNA was proposed 40 years ago [Q . 
In recent years its conducting properties were tested in 
several contexts: DNA has been used as a template on 
top of which nanowires are built j2j and there is hope that 
self- assembled DNA junctions between conductors will 
function soon. There has also been renewed interest in 
building and characterizing molecular wires of plain DNA 
between two metal leads j^, ^, |], ^, ^ , possibly combined 
with changes of the major constituents of the molecule 
to improve its properties . There is ongoing debate as 
to whether DNA is conducting, semi-conducting, insu- 
lating or all of the above p. There are several reasons 
for this controversy. There is great diversity of the DNA 
forms in terms of its composition, length, and structure. 
Counterions and impurities that can attach to the phos- 
phate groups and the grooves of the molecule need to be 
taken into account. Finally, environmental factors such 
as temperature and contact resistance can play a major 
role. 

Experiments done 50 years ago by Wilkins et al. [p^ , 
suggested that overstretched DNA (that is, substan- 
tially longer than its natural length) undergoes a tran- 
sition to a structure that can accommodate elongation 
up to twice the length of relaxed DNA. A major break- 
through in understanding the mechanical properties of 
DNA was achieved through single molecule stretching ex- 
periments In, 12, 13|. It was shown that by gradually in- 
creasing the axial stretching force F along the DNA helix, 
the molecule first uncoils, then, at 5 pN, it reaches its 
natural length, and then exhibits a stiff elastic response 
for forces of up to 50 pN. At F 50 to 80 pN, DNA under- 
goes a pronounced and abrupt structural transformation 
to a yet unknown structure that is elongated by more 
than 50%. Further increase in force leads to breaking 
after it has stretched to twice its natural length. In this 
letter we argue that these overstretched forms of DNA 
are relevant in the quest for novel device components, 
we develop a visualization method for understanding the 



nature of electronic states in these forms, and we present 
the first detailed discussion of their electronic properties. 

We use an efficient quantum mechanical electronic 
structure method |lj] that treats all the valence electrons 
and that accounts for the charge transfer involved in 
biological molecules and semiconducting surfaces. This 
approach, called self-consistent-charge density-functional 
tight-binding method, is based on an expansion of the 
Kohn-Sham energy functional JTsI , in terms of the 
charge density. The zeroth order term results in a tight- 
binding Hamiltonian, and the second order term incor- 
porates the charge redistribution. The method is able to 
deal with the very large numbers of atoms involved in bi- 
ologically relevant systems. As a test, we performed cal- 
culations on a periodic homopolymer of poly(C)-poly(G) 
with 10 base pairs per turn and reproduced the results of 
dePablo et al. Q , which to our knowledge is the only ab 
initio electronic structure calculation of a periodic DNA 
structure [Q. Our results are in excellent agreement to 
their calculations: we find a 2.1 eV bandgap between 
the highest occupied and lowest unoccupied eigenstates, 
compared to 2.0 eV reported in Ref. The wavefunc- 
tions of the highest occupied valence band are located 
on the guanines while the wavefunctions of the lowest 
unoccupied band are located on the cytosines, also in 
agreement with Ref. ||]. 

For the purpose of the present study, calculations 
were performed on a variety of DNA polymers with 
different compositions employing both infinite periodic 
and finite structures. We limit our discussion to 
poly (CG) -poly (CG) as a representative structure for 
which well established overstretched configurations have 
been published jlSj . Our results on finite DNA polymers 
show the existence of states localized at the opposite ends 
of the molecule. The energy levels of these states lie be- 
tween the regular bands that appear in the periodic struc- 
tures. We removed these states before analyzing the elec- 
tronic properties, and will refer to the remaining states 
as the "restricted band structure" of the non-periodic 
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configuration. As a test, we compared 18 base pair re- 
stricted band sturcture with periodic structures with 10 
base pairs per unit cell and found the band structures 
virtually indistinguishable. A boundary region of 3 base- 
pair layers at each end of the molecule was sufficient to 
contain the spurious end states. Fig. |^ displays the band 
structure of a periodic sequence with a unit cell consist- 
ing of 10 base pairs. The sharp features of the density of 
states are due to the one dimensional nature of DNA. In 
all periodic sequences there is an extra symmetry, which 
for the poly(CG)-poly(CG) structure is the translation 
by two base pairs along the stack, together with a rota- 
tion of 72 degrees around the helical axis: this symmetry 
is generated by the screw operator that is isomorphic to a 
pure translation operator. Structures obtained by short 
(50 ps) molecular dynamics runs at 50 K using classi- 
cal force fields ||l^, indicate that minor deformations of 
the ideal configuration destroy the sharp features in the 
electronic structure. 

We now move on to a discussion of overstretched DNA. 
First we elaborate on how this is relevant in the con- 
struction of DNA nanowires. One way to attach DNA 
molecules to substrates is the technique of "molecular 
combing" Q. This method uses the meniscus forces 
that develop between the solution surface and the device 
surface template to uncoil the DNA and stretch it be- 
tween two electrodes as the device is pulled out of the 
solution. A description of recent experiments estimates 
that the forces that are exerted on the DNA molecules 
are of order 100 pN |21 , more than enough to stretch 
them past their natural length. This is indicated by re- 
cently built devices |^ where the distance between the 
electrodes was about 30% longer than the natural length 
of the DNA molecules. The overstretching can have a 
variety of consequences on the electronic states. 

Elongated DNA structures have been determined 
in the pioneering study of Lebrun and Lavery [ p^ , 
which model the adiabatic elongation of selected DNA 
molecules. There are two distinct modes of overstretch- 
ing corresponding to pulling the opposite 3'-3' ends and 
5'-5' ends. The former mode leads to DNA unwinding, 
consistent with a recent estimate of the helicity of over- 
stretched DNA whereas the latter leads to a con- 
traction of the diameter of the helix. Both overstretch- 
ing modes can accommodate elongations of up to 90% 
without breaking the molecule. 

We next present the inherent electronic states of 
the overstretched structures and then consider electron 
transport aspects. For all the calculations we used acidic, 
dry, 18 base-pair strings of DNA, and we performed the 
projection of the end states in the manner stated before. 
The quasi one-dimensional structure of all the forms leads 
to the formation of clear minibands in the energy spec- 
trum, similar to those shown in Fig. |l|. The band struc- 
ture of the overstretched forms have an overall resem- 
blance to the band structure of the original unstretched 
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FIG. 1; The band structure around the Fermi energy (left) 
with the corresponding density of states (right) for the peri- 
odic poly (CG) -poly (CG) form. The inset shows a detail of 
the density of states at the highest occupied valence band. 
We denote the bandgap around the Fermi energy (g^), the 
widths of the highest valence (w") and lowest conduction (wj) 
minibands, and the gaps between these and the neighboring 
minibands (gi and gi). 



form, and many of the minibands near the Fermi en- 
ergy can be continuously mapped to their counterpart 
in the unstretched configuration. This is mostly because 
of the molecular structure of the overstretched forms as 
will become clear in the following discussion. For some of 
the most deformed structures, especially the 90% over- 
stretching in the 5'-5' mode, several of the minibands are 
mixed and small but clear gaps present in less deformed 
structures disappear. 




FIG. 2: The left panel shows a piece of poly (CG) -poly (CG) 
B-DNA. Along its helical axis is a cylinder on which we plot 
the contours of the valence band electron density. The right 
panel shows the same region for a 30 % stretched structure. 

We have devised a new scheme for visualizing the elec- 
tronic states of DNA, which naturally accounts for the 
symmetry of this molecule. The basic concept of this 
scheme is to unfold contour plots of the band electron 
density evaluated at cylindrical surfaces centered on the 
helical axis, as shown in Fig. ^. In this figure we show the 
structure at its natural length together with one that is 
elongated by 30% in the 3'-3' overstretching mode. With 
the help of these plots we can understand most of the 
band structure features of the stretched forms. The elon- 
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gation to the overstretched form is achieved by changing 
the dihedral angle configuration of the backbone of the 
molecule, thus leaving the local part of the orbitals es- 
sentially intact. 

In Fig. ^ (a) we present the valence band electron den- 
sity of selected stretched poly (CG) -poly (CG) structures 
using our visualization scheme. We choose the radii of 
the cylinders to follow the atom with the largest electron 
density for the band that we present. In the unstretched 
structure, the 7r-band is along the guanines. Upon elon- 
gation, orbitals from the picture frame are removed while 
the remaining ones rotate partially as dictated by the new 
orientation of base pairs. The poly(CG)-poly(CG) struc- 
ture contains pairs of guanines that are closer to each 
other. These guanines display an overlap of the elec- 
tronic densities as seen in the top left frame of Fig. |^ (b) . 
Elongation in the 3'-3' mode keeps such paired guanines 
close to each other which is reflected in the overlap in 
their electronic densities; the overlap of the electronic 
densities between guanines from different pairs decreases 
with elongation. In Fig. ^ (b) we show the conduction 
band electron density. Note again how the orbitals ro- 
tate as the structure is being overstretched, following the 
rotation of the bases. Note that in the extreme 90% 
stretching case, the orbitals become perpendicular to the 
helical axis. 

The changes in the structural features should be of 
profound consequence to any model dealing with elec- 
tron transfer since, to a first approximation, the over- 
lap matrix elements depend exponentially on the dis- 
tance between orbitals. Because of the tt character of 
the states near the Fermi level, a non-parallel alignment 
will strongly affect the values of the overlap integrals. In 
order to quantify these effects, we considered the nearly 
degenerate eigenstates in the first valence miniband. We 
used the linear combination of these states that gives nor- 
malized, orthogonal states which are maximally lo- 
calized along the helical axis within the subspace spanned 
by the delocalized states. The resulting valence states are 
localized on guanine bases. With these states, we calcu- 
lated the hopping matrix elements: 

t,, = -i(vl/,|V2|vI/^-), (1) 

between nearest-neighbor sites. Note that for the 
poly (CG) -poly (CG) structure every guanine has 2 gua- 
nine neighbors, one above and one below as indicated in 
Fig. ^, separated by distances of 5.3 A and 4.2 A. The 
larger of the two matrix elements corresponds to hopping 
between sites that are close together; its value remains 
roughly constant for the 30% overstretched structures. 
The smaller of the two matrix elements, which corre- 
sponds to sites that are further apart, is dramatically 
reduced for all the overstretched forms. Since the two 
different hopping terms are connected in a 1-D series the 
smallest term will determine the bottleneck for electron 



transport. The absolute square values of these small hop- 
ping matrix elements between these neighbors are given 
in Fig. ^ (a). We infer from this figure that if the con- 
duction mechanism is through holes there should be a 
dramatic drop in conductivity with 30% overstretching, 
which in the 3'-3' mode drops further with elongation. In 
Fig. 1^ (b) we give the hopping matrix elements between 
states in the conduction band. In the cases where the 
pairing between nearest neighbor pairs have very differ- 
ent matrix elements we only list the one that would be 
the bottleneck for this channel. Due to the stronger local- 
ization of the conduction band and because of the rather 
large distance between the occupied sites, the hopping 
matrix elements are substantially smaller than the ones 
between the valence states in the structure of natural 
length. In fact, it is the site that is two bases away which 
has the strongest overlap in this structure: the C6 atoms 
between these bases are actually closer to each other (9.0 
A) than the C6 atoms between consecutive base pairs 
(10.0 versus 10.3 A). The first and second neighbor hop- 
ping terms are indicated in Fig. || (b). Since these terms 
are connected in parallel, this time it is the maximum 
of the two that will play the dominant role in the con- 
duction mechanism and it is this one that is included in 
Fig. H (b). We note that if conduction takes place by 
electrons the conductivity drops with overstretching in 
the 3'-3' mode but not necessarily in the 5'-5' mode. 

The concept of overstretched DNA in relation to elec- 
tronic device components appears relevant to a variety 
of experimental and theoretical issues. Experimentally 
it is possible, and sometimes even unavoidable, to create 
overstretched DNA molecules. The window of allowable 
forces that can be applied to stretch DNA to its natural 
length, but not overstretch, is small (from 5 pN to 50 pN) 
compared to the forces applied during molecular comb- 
ing. Another important aspect is owed to the nature of 
the extreme stretched forms, it seems possible to expose 
the bases to the exterior of the molecule allowing impuri- 
ties to enter the crucial base-pair stack region; this could 
result in controlling the electronic properties with doping. 
This might also enable direct visualization of parts of the 
electronic wavefunctions which was is not possible in the 
natural DNA forms because of the insulating backbone 
layer. On the theoretical level, the major existing mod- 
els (superexchange [Q, hopping and polaron psf), 
of pure DNA conductivity will be drastically affected by 
the consequences of overstretching, including changes in 
the hopping matrix elements and in the low-energy vi- 
brational modes. For long fibers of homogeneous DNA 
(either a homopolymer, or a short repeated sequence) 
we expect stretching to appear homogeneously. Regions 
of varying richness in AT or GC mixture would induce 
different lengthenings at different parts of the molecule, 
permitting conduction mechanisms to change along the 
molecule. 

We thank Leo Kouwenhoven and Nina Markovic for 
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FIG. 3: Cylindrical contour map of the highest valence (a) and the lowest conduction (b) band electron density, for the 
poly(CG)-poly(CG) DNA (top left) and the 30%, 60%, and 90% overstretched forms (right). The bottom left panel of each 
figure shows a drawing of the CG base pairs along the cylinder; the arrows indicate bases with the largest couplings. In Fig. (a) 
the 90 % overstretched 5'-5' form is not shown because the highest valence band mixes with the lower bands. The vertical 
axis covers a fixed region corresponding to 10 base pairs in the unstretched form. The horizontal axis runs from to 27r. The 
radius of the cylinder follows the guanine C5 atom in (a), and the cytosine C6 in Fig. (b). The cylinder surface unit element is 
kept fixed: this results in a variation of the length of the horizontal axis, proportional to the cylinder radius. The color coding 
(shown in the color bar between the unstretched and stretched graphs) is fixed between different structures and covers 5 orders 
of magnitude. 
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FIG. 4: The values of \tijf from Eq. (1) for (a) the valence 
and (b) the conduction band. The values are given eV^ mul- 
tiplied by a factor of lO'^^ to facilitate relative comparisons, 
for overstretching of 0, 30, 60, 90 % relative to the natural 
length. The blue (left) column at each entry corresponds to 
the 3'-3' stretching mode and the red (right) column to the 
5'-5' mode. In Fig. (a), the entry at the 90 % stretching of 
the 5'-5' mode is missing because of the mixing the valence 
band with lower bands. 



sharing their experimental insight with us. We are grate- 
ful to Richard Lavery for providing the overstretched 
structures. We thank Hatem Mehrez for his critical read- 
ing of the manuscript. This work was supported in part 
by Harvard's Materials Research Science and Engineering 
Center, which is funded by the National Science Foun- 
dation. This work was partially supported by National 
Computational Science Alliance under MCB010005N and 
utilized the Boston University Scientific Computing and 
Visualization facilities. RLB was supported by a national 
science foundation graduate research fellowship. 



[1] D. D. Eley and D. I. Spivey, Trans. Faraday Soc. 58, 411 
(1961). 

[2] E. Braun, Y. Eichen, U. Sivan, and G. Ben-Yoseph, Na- 
ture 391, 775 (1998). 



[3: 
[4: 

[5] 
[6] 
[7 

[s: 

[9 
[10 

[11 

[12 
[13 



[14 

[is: 
[16: 

[17 



H.-W. Fink and C. Schonenberger, Nature 398, 407 

(1999) . 

D. Porath, A. Bezryadin, S. De Vries, and C. Dekker, 
Nature 403, 635 (2000). 

P. de Pablo et al, Phys. Rev. Lett. 85, 4992 (2000). 
A. Kasumov et al, Science 291, 280 (2001). 
K.-H. Yoo et al, Phys. Rev. Lett. 87, 198102 (2001). 
A. Rakitin et al, Phys. Rev. Lett. 86, 3670 (2001). 
H. W. Fink, Cell. Mol. Life Sci. 58, 1 (2001). 
M. H. F. Wilkins, R. G. Gosling, and W. E. Seeds, Nature 
167, 759 (1951). 

S. B. Smith, Y. J. Cui, and C. Bustamante, Science 271, 
795 (1996). 

P. Cluzel et al. Science 271, 792 (1996). 
T. R. Strick, J. F. AUemand, D. Bensimon, and V. Cro- 
quette, Annu. Rev. Biophys. Biomolec. Struct. 29, 523 

(2000) . 

M. Elstner et al, Phys. Rev. B 58, 7260 (1998). 
P. Hohenberg and W. Kohn, Phys. Rev. B 136, B864 
(1964). 

W. Kohn and L. J. Sham, Phys. Rev. 140, 1133 (1965). 
There exist early papers on computational modeling of 
simplified DNA structures or of small, finite piecies of 
DNA. Some early work of Ladik and collaborators, is pre- 
sented in A. K. Bakhshi, P. Otto, J. Ladik, and M. Seel, 
Chem. Phys. 108, 215, (1986), and references therein. 
[18] A. Lebrun and R. Lavery, Nucleic Acids Research 24, 
2260 (1996). 

[19] W. D. Cornell et al, J. Am. Chem. Soc. 117, 5179 (1995). 
[20] A. Bensimon et al, Science 265, 2096 (1994). 
[21] D. C. G. Klein et al, Appl. Phys. Lett. 78, 2396 (2001). 
[22] M. W. Wu, L. P. Kouwenhoven, and L. L. Sohn (unpub- 
lished) . 

J. F. Leger et al, Phys. Rev. Lett. 83, 1066 (1999). 
J. Jortner, M. Bixon, T. Langenbacher, and M. E. 
Michel-Beyerle, Proc. Natl. Acad. Sci. 95, 12759 (1998). 
P. T. Henderson et al, Proc. Natl. Acad. Sci. U. S. A. 
96, 8353 (1999). 



[23] 
[24] 

[25] 



